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Abstract — We consider a MIMO fading broadcast channel and 
compare the achievable ergodic rates when the channel state 
information at the transmitter is provided by "analog" noisy 
feedback or by quantized (digital) feedback. The superiority 
of digital feedback is shown, with perfect or imperfect CSIR, 
whenever the number of feedback channel uses per channel 
coefficient is larger than 1. Also, we show that by proper design 
of the digital feedback link, errors in the feedback have a minor 
effect even by using very simple uncoded modulation. Finally, 
we show that analog feedback achieves a fraction 1 2F of 
the optimal multiplexing gain even in the presence of a feedback 
delay, when the fading belongs to the class of "Doppler processes" 
with normalized maximum Doppler frequency shift < F < 1/2. 

I. Model setup and background 

We consider a multi-input multi-output (MIMO) Gaussian 
broadcast channel modeling the downlink of a system where 
the base station (transmitter) has AI antennas and K user 
terminals (receivers) have one antenna each. A channel use 
of such channel is described by 

yk^h."x + Zk, k = l,...,K (1) 

where yk is the channel output at receiver k, Zk ^ C3\f(0, A^o) 
is the corresponding AWGN, € C*^ is the vector of 
channel coefficients from the fc-th receiver to the transmitter 
antenna array and x is the channel input vector. The channel 
input is subject to the average power constraint E[|xp] < P. 

We assume that the channel state, given by the collection 
of all channel vectors H = [hi, ... , h/<-] G C^^^^, varies in 
time according to a block fading model where H is constant 
over each frame of length T channel uses, and evolves from 
frame to frame according to an ergodic stationary jointly 
Gaussian process; i.i.d. block-fading channel, where the entries 
of H are Gaussian i.i.d. with elements ~ CINf(0, 1) is a special 
case of this. 

A. Capacity results 

If H is perfectly and instantaneously known to all terminals 
(perfect CSIT and CSIR), the capacity region of the channel 
([T]i is obtained by MMSE-DFE beamforming and Gaussian 
dirty-paper coding (see [1,2] and references therein). Because 
of simplicity and robustness to non-perfect CSIT, simpler 
linear precoding schemes with standard Gaussian coding have 
been extensively considered. A particularly simple scheme 
consists of zero-forcing (ZF) beamforming, where the transmit 



signal is formed as x = Vu, such that V G C ^ is a zero- 
forcing beamforming matrix and u e C^^ contains the symbols 
from K independently generated Gaussian codewords. For 
K < M, the k-th column Vfc of V is chosen to be a unit 
vector orthogonal to the subspace §k = spanjhj : j ^ k}. In 
this case, the achievable sum rate is given by 

.na. Ve [log fl+W-'f 

^^E[VUH)]<P ^ L V No J 

(2) 

We consider the situation where K — M, and thus do not con- 
sider user selection. Furthermore, we are mainly interested in 
the high-spectral efficiency regime, where we can characterize 
the achievable sum rate as KlogP/No + 0(1), and k is the 
"system multiplexing gain" or "pre-log factor" of the ergodic 
sum rate. Hence, it is well-known that using uniform power 
Vk = P/M for all k — 1,...,M, rather than performing 
optimal water-filling, incurs a loss only in the 0{\) term, and 
we shall restrict to this choice in the rest of this paper 

It is well-known that, under perfect CSIT and CSIR, both 
the optimal "Dirty-Paper" sum-rate C and the zero-forcing 
sum-rate R^^ are equal to M log P/Nq + 0(1). On the 
contrary, under non-perfect CSIT the rate sum may behave 
in a radically different way; for example, if there is perfect 
CSIR and no CSIT when H has i.i.d. Gaussian entries, the 
sum rate is equal to XogP/N^ + 0(1) [1] 

B. Channel state feedback models 

We consider some specific CSIT and CSIR models and 
derive lower-bounds to the corresponding achievable ergodic 
rates by analyzing a naive beamforming scheme that computes 
a mismatched ZF beamforming matrix V from the CSIT. In 
particular, we consider an "analog" CSIT feedback scheme 
where the transmitter observation at frame time t is given by 

{G(r) = V^H(t) + W(r) : r = -c5o, . . . ,t ~ d] (3) 

where {W(r)} is a spatially and spectrally white Gaussian 
process with elements ~ C3Nf(0, N^) and d is the feedback 
delay. This models the case where the channel coefficients are 
explicitly transmitted on the reverse link (uplink) using un- 
quantized quadrature-amplitude modulation [3-6]. The power 
scaling (3 corresponds to the number of channel uses per 
channel coefficient, assuming that transmission in the feedback 
channel has fixed peak power P and that the channel state 
vector is modulated by a (3MxM unitary spreading matrix [3]. 



A simplifying assumption of this work is that we consider no 
fading and orthogonal access in the CSIT feedback link, and 
we assume that the SNR on the feedback channel is equivalent 
to the un-faded downlink SNR (P/Nq). 

A different CSIT feedback approach is based on quantizing 
the channel vector at each receiver and transmitting back to the 
base station a packet of B bits, representing the corresponding 
quantization index. If a random ensemble of quantization 
schemes is used (referred to as Random Vector Quantization, 
or RVQ), in [7, Theorem 1] it is shown that the gap between 
ZF with ideal CSI and the naive ZF scheme is given by 



Ai?quant. <log(l + ^2- 



(4) 



II. Rate gap bound for analog CSIT feedback 

In the case of i.i.d. block fading and no feedback delay, 
the analog CSIT feedback yields the observation of G = 
\//3PH + W at the beginning of every frame. The transmitter 
computes the MMSE estimate of the channel matrix, H = 
Na+pp ^' ^'^^ column Vfc of V is a unit vector orthogonal 
to the subspace Sfc = spanjhj : j ^ k}. Notice that we can 
write H = H + E, where H and E are mutually independent 
and have Gaussian i.i.d. components with mean zero and 
variance and a'^ — (1 + PP/Nq)^^, respectively. 



The signal at the A:-th receiver is given by 



Vk 



Zk 



(5) 



We assume that the frame duration is long enough such 
that some training scheme can be used in the downlink 
channel. Training allows each receiver to estimate: 1) the 
useful signal coefficient, — (h^Vfe) and 2) the variance 
of the interference plus noise (k = + -^fc' 



given by Sfe = E |CfcP|efc,H = TVq + ^ 



P/M. 



This conditioning is due to the fact that T,k is estimated on 
each frame, and the coefficients (e^Vj) are constant over 
each frame and change from frame to frame, following the 
block i.i.d. fading model. The maximum achievable rate of 
user k subject to the above assumptions is lowerbounded by 
assuming a Gaussian input Uk ^ u'j^ ^ G'N{Q, P/M), and by 
considering the worst-case noise plus interference distribution 
in every frame. Using stationarity and ergodicity, we have0 



Rk 



> 



E 



E 



log 1 



(6) 



where (a) follows from [8], noticing that a^uf and are 

uncorrected (even after conditioning on ak, Sa). 

Next, we shall bound the rate gap incurred by the naive 
ZF beamforming and analog feedback with respect to the ZF 

'With some abuse of notation, tlie term in the second line of (6) have the 
following meaning: 

^ [^°fCfc:E[|CfcPl<Efc /(n^;j/fc|afc,Sfc)] 

I ^^^(:k-n\ikF]<<^ ^^^k'yk\ak,^k = a-)dF{(T) 
where F{cr) denotes the cdf of S^. 



beamforming with ideal CSIT. Denoting by Rf^ the rate of 
user k with uniform (across users) and constant (in time) 
power allocation ^^(H) = P/M in ^, we have 



'analog 



Rk 



log 
log 



1 + 



1 + 



-E 



log 1 + 



-E 



+E 



log 1 + 



NoM 
NqM 

(E,^JeHv,|2 + |a,|2^ 



NqM 



log 1 



E 



(a) 
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(b) 
< 



1 + E 



NqM 



log 1 + 



NoM 



log 1 + 



o-2p M - 1 



(7) 



where (a) follows from the fact that J2j^k l^fe^jP + lofcP 
stochastically dominates jhjjvfcp since |afcp and jh^Vfep are 
identically distributed, (b) follows from Jensen's inequality 
and the final expression (c) follows by noticing that the V is a 
deterministic function of H and therefore it is independent of 
E. Therefore, we can write E[|eJ^Vjp] = E[v^E[eke\^]vj] = 
(TgE[|vjp] = (Tg, since Vj has unit norm by construction. 

III. Comparison with quantized CSIT feedback 

In this section we compare analog and digital feedback 
under the assumptions of perfect CSIR, no feedback errors, 
and no feedback delay. Replacing the estimation error variance 

a2 = (1 + pP/No) 
obtain: 



in (Q and further upper bounding we 



AR 



analog 



(8) 



Let us now consider digital feedback over the same channel. 
The rate gap obtained in [7, Theorem 1] and reported in (HI is 
further upperbounded by log(H-(P/A^o)'2^~). Let us assume 
(very unrealistically) that the digital feedback link can operate 
error-free and at capacity, i.e., it can reliably transmit log(l + 
P/Nq) bits per symbol. For the same number of feedback 
channel periods, f3M, the number of feedback bits per mobile 
is B = /3Mlog2(l + P/Nq). Replacing this into the rate gap 
bound, we obtain: 



AR. 



quant, 



< log 1 



P/No 



(9) 



(1 + p/Nor 

If /? = 1 the quantized and analog feedback achieve essentially 
the same rate gap of at most 1 b/s/Hz. However, if /3 > 1, 
unlike the analog feedback case, the rate gap of the quantized 
feedback vanishes for P/Nq oo. and digital is far superior 
to analog for (3 > 1. 

This conclusion finds an appealing interpretation in the 
context of rate-distortion theory. It is well-known (see [9] 
and references therein) that analog transmission is an optimal 
strategy to send a Gaussian source over a Gaussian channel 



with minimal end-to-end quadratic distortion. In our case, the 
source is the Gaussian channel vector and the noisy channel 
is the feedback AWGN channel with SNR P/Nq. Hence, the 
fact that analog feedback cannot be essentially outperformed 
for /3 = 1 is expected. However, it is also well-known that 
if the channel rate is larger than the source rate (i.e., less 
than one Gaussian source symbol arrives per channel symbol, 
which corresponds to /3 > 1 in our case), then analog is 
strictly suboptimal as compared to separate source and channel 
coding because the distortion with analog transmission scales 
as 1//3 whereas it decreases exponentially with /3 (i.e., along 
the vector quantizer R-D curve) for digital transmission. 

IV. Effects of Imperfect CSIR 

We now consider the scenario where each receiver has only 
a noisy estimate of its channel acquired via downlink training. 
In order to allow for channel estimation, piM shared pilots 
(/?! > 1 symbols per antenna) are transmitted. Each receiver 
estimates its channel on the basis of Y = ^//3iPH + Z, which 
yields (after MMSE estimation) Gaussian error with variance 
(1 + PiP/Nq)^^ . Terminals feed back channel information 
immediately after completion of this training phase. After the 
transmitter has chosen beamforming vectors on the basis of the 
channel feedback, an additional round of downUnk training 
is performed to enable coherent detection and allow each 
terminal to estimate its useful signal coefficient ak — h^Vk. 
This can be accomplished in f32M symbols by transmitting 
along each of the beamforming vectors for (32 symbols. If 
MMSE estimation of is performed, we have ak — ak + fk 
where fk and cik are independent complex Gaussian's with 
variance aj.- — jy^jj^J^^p and 1 — erf, respectively]! 

Under this set of assumptions, a lower bound to 
I{uk', ykl^k) can be derived using techniques similar to those 
in [10, 11]. Using this lower bound and some steps similar to 
those leading to (|7]i, the following upper bound to the rate gap 
can be reached at: 

Ai? < log2(l + ^(a2. + (M-l)EKv,f]),^ 

where the multi-user interference term E[|h^Vjp] depends on 
the CSIT and thus on the channel feedback (/3) as well as 
the accuracy of the initial training (/3i). We again assume that 
f3M symbols are devoted to channel feedback (per mobile). If 
analog feedback is used, we get an upper bound of: 

Ai?.„.o. < log(l + ^ + ^ + i) (10) 

In the case of digital feedback, under the assumption that B = 
/3Mlog2(l + P/Nq) feedback bits per mobile are sent in an 
error-free manner, we get: 

Comparing ( fTOl i and ( fTTT i we come to the same general 
conclusions as in Section |III1 if /3 = 1 then digital and analog 
are equivalent, but if /? > 1 digital is superior to analog 
because the effect of feedback noise vanishes at high SNR 
for digital but does not do so for analog. 

-Note that additional training is required because terminals do not know the 
channels of other terminals, and thus are not aware of the chosen beamforming 
vectors. 
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Fig. 1. Quantized vs. Analog Feedback with Imperfect CSIR. 



There are, however, some important differences with the 
perfect CSIR scenario. First note that the imperfect CSIR leads 
to residual interference that does not vanish with SNR; as 
a result, the rate gap is not driven to even when f3 > 1, 
assuming /3i and P2 are fixed. In addition, when /3i « /32 ~ 1, 
imperfect CSIR seems to have a considerably stronger effect 
than feedback noise, thereby reducing the magnitude of digital 
feedback's advantage. These effects are both visible in Fig. 
[T] where analog and digital feedback curves are plotted for 
Pi = p2 = ^ and (3=1 and (3 = 2, along with the throughput 
of an imperfect CSIR/perfect FB system. 

Finally we comment on the tradeoff between downlink 
training (/3i) and channel feedback {(]). Since downlink pilots 
are shared, training consumes only (3iM channel symbols. 
Channel feedback, on the other hand, requires (3M channel 
symbols per mobile. If the M terminals can simultaneously 
transmit on the feedback channel, perhaps utilizing the M- 
antenna receive array at the base as described in [3], then (3 
and (3i are equivalent in terms of system resources. For the 
case of analog feedback, from ( fTOl i we see that (3 and (3i 
should be chosen equal. On the other hand, if digital feedback 
is used, it is only necessary to choose (3 > 1 (so that the 
effect of feedback noise vanishes), and the remainder of the 
resources should be devoted to downlink training, i.e., to (3i. 
This is an additional advantage to digital whenever (3+(3i > 2. 
Note that there is also a tradeoff between /3i and (32, but that 
the effect of the initial training ((3i) is considerably stronger 
than the second phase. 

V. Effects of CSIT feedback errors 

We now investigate the impact of removing the optimistic 
assumption that the quantized feedback channel can oper- 
ate error-free at capacity. We consider a very simple CSIT 
feedback scheme that certainly represents a lower bound 
on the best quantized feedback strategy. The user terminals 
perform quantization using RVQ and transmit the feedback 
bits using simple uncoded QAM. No intelligent mapping of the 
quantization bits onto the QAM symbols is used, and therefore 
even a single erroneous feedback bit from user k results in 
CSIT that is completely independent (due to the properties 
of RVQ) of the actual fc-th channel vector Since uncoded 
QAM is used, error detection is not possible and the base 
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Fig. 2. Quantized feedbacl: with QAM modulation. 



Station computes beamforming vectors based on the possibly 
erroneous feedback. 

We again use fiM symbol periods to transmit the feed- 
back bits. There is a non-trivial tradeoff between quanti- 
zation and channel errors. In order to maintain a bounded 
gap, feedback must be scaled at least as [M — l)log2(l + 
P/Nq) « AI \og2 P/Nq. Therefore, we consider sending B ~ 
aM log2 P/No for 1 < a < /? bits in /3Af symbol periods, 
which corresponds to ^ log2(i-'/A*'o) bits per QAM symbol. 

From [12], using the fact that the QAM constellation size is 
equal to L = (P/Nq)^ , we have the following upper bound 
to the symbol error probability for QAM modulation: 



1-Q//3N 



(12) 



For a — (3 (which means trying to signal at capacity with 
uncoded modulation!) Pg does not decreases with SNR and 
the system performance is very poor However, for a/ (3 < 1, 
which corresponds to transmitting at a constant fraction of 
capacity, Ps ^ as P/Ng oo. The upper bound on the 
error probability of the whole quantized vector (transmitted 
in (3M symbols) is given by Pejb = 1 - (1 - Ps)'^'^' ■ A 
lower bound on the achievable ergodic rate is obtained by 
assuming that when a feedback error occurs for user k its 
SINR is zero while if no feedback error occurs its rate is given 
— Ai?quant., that is, the rate of ideal ZF decreased by the 
(upper bound to) the rate gap. It follows that the ergodic rate 
of user k is upperbounded by 



VI. Effects of CSIT feedback delay 

We consider now the case of analog feedback (assuming 
perfect CSIR) when each entry of H evolves independently (in 
the block-fading way described earlier) according to the same 
complex circularly symmetric Gaussian stationary ergodic 
random process, denoted by {h{t)}, with mean zero, variance 
1 and power spectral density (Doppler spectrum) denoted by 

5,(0, [-1/2,1/2]. 

Because of stationarity, without loss of generality we can 
focus on t = 0. We are interested in the linear MMSE 
estimation of h{t) from the observation {g{T) : r = —oo,t — 
d} where, following the analog feedback model (|3]l, we let 
g{T) = h{T) + w{T), with w{t) i.i.d. - eJ<{0,S) and S ^ 0. 
In particular, we consider the case of 1-step prediction (d = 1) 
and the case of filtering (d = 0). From classical Wiener 
filtering theory [13], we have that the prediction error is given 
by 

eiiS) = exp y'J^ \og{5 + 5„(0)dej - S (14) 
and that the filtering MMSE is given by 

eo(<5) - 



(15) 



We shall discuss the rate gap bound (O letting al — 
^d{N{)/ {I3P)) for d = 0,1, under different assumptions on 
the fading process {h{ty}. We distinguish two cases: Doppler 
process and regular process. We say that is a Doppler 

process if Sh{^) is strictly band-limited in [—F,F], where 
F < 1/2 is the maximum Doppler frequency shift, given 
by F = ^Tf, where v is the mobile terminal speed 
(m/s), fc is the carrier frequency (Hz), c is light speed (m/s) 
and Tf is the frame duration (s). Furthermore, a Doppler 
process must satisfy J^^log Sh{£,)d^ > ~oo. Following [14], 
we say that {h{t)} is a regular process if ei(0) > 0. In 
particular, a process satisfying the Paley- Wiener condition [13] 

I-i%^ogSh{Od^ > -oo is regular. 

A Doppler process satisfying our assumptions has prediction 



1-2F . 



cxp 



log{S + ShiO)dn - 6 (16) 



No feedback delay (d 



Pai 



0). In this case 
iVo^(t) 



(17) 



Rk > (1 - Psf''' {RT - log (1 + (^/^o)'"")) (13) Hence, limp^oo P<t, 



No 



Choosing 1 < a < (3 we achieve both vanishing Pg and van- 
ishing Ai?quant. as P / Nq OO. Thus, cvcn under this very 
simple CSIT feedback scheme the optimal ZF performance 
can be eventually approached for sufficiently high SNR. 

Fig. |2] shows the ergodic rate achieved by ZF beamforming 
with quantized CSIT and QAM feedback transmission for 
M = K = A, independent Rayleigh fading, /3 = 4 and 
different values of a. It is noticed that by proper design of 
the feedback parameters the performance can be made very 
close to the ideal CSIT case. 



^ for both Doppler and regular 
processes. For the latter, this is clear from the fact that 
ei(0) > 0. For the former, this follows from (fTST l. Applying 
Jensen's inequality and the fact that / Sh{(,)d£^ — 1, we arrive 
at the upper bound 



1 

2F 



m 



(18) 



Using the fact that log is increasing, we arrive at the lower 
bound 



cxp 



(19) 




These bounds yield that ei(A^o//3-P) = nP-'^^-'^^'^ +0{ll P) 
for some constant k. Hence, ei = ©(P^*^^^^^^) while 5 = 
0{l/P), and the hmits holds. 

We conclude that in the case of no feedback delay the 
estimation error is essentially dominated by the instantaneous 
observation and not much improvement can be expected by 
taking into account the channel memory if analog feedback 
is used. With quantized feedback the same may not be true 
because it is possible to exploit memory by feeding back only 
the innovation process [15]; this is under investigation. 

Feedback delay {d = 1). In this case, the behavior of 
Doppler versus regular processes is radically different. For 
Doppler processes, using (fTsl l and (fT9] l, we have that Pa^ = 
Pei{Na/l3P) = kP'^^ + 0(1). It follows that the achievable 
rate sum is lowerbounded by 



M 

E 

k=l 



Rk > M{1 - 2F)logP + 0(l) 



(20) 



which implies a multiplexing gain of M{1 — 2F). 

For regular processes, on the contrary, we have that Pa^ > 
Pei(O) = 0{P). Hence, the rate gap grows like logP and 
the achieved multiplexing gain is zero. Furthermore, it can be 
shown that the following is actually an upper bound to the 
per-user rate, even when the feedback is noiseless: 

1 \ ^(M) 



Rk < log2 



1 



1 



+ (Af - 1) 



loge2 



1 



1 



logg2 V2Af-l ' 2Af-2/ "-^^^ 

In conclusions, the most noteworthy result of this analysis 
is that under common fading models (Doppler processes), the 
analog feedback scheme achieves a potentially high multi- 
plexing gain even with realistic, noisy and delayed feedback. 
Notice for example that with mobile speed ?; = 50 km/h, 
/c = 2 GHz, and frame duration 1 ms, we have F — 0.0926. 
With M = 4 antennas we achieve a yet respectable pre-log 
factor equal to 3.26 instead of 40 

^It is interesting to notice here the parallel with the results of [14] on the 
high-SNR capacity of the single-user scalar ergodic stationary fading channel 
with no CSIR and no CSIT, where it is shown that for a class of non-regular 
processes that includes the Doppler processes defined here, the high-SNR 
capacity grows like log P, where L is the Lebesgue measure of the set 
{5 e [-1/2, 1/2] : ShiO = 0}. In our case, it is clear that li = l-2F. 



Figs. [3] and |4] show the achievable ergodic rates for the 
Jakes' "Jo" correlation (strictly band-limited) and the Gauss- 
Markov AR-1 correlation (regular process) for different first- 
lag correlation values. For the AR-1 process with d = 1 
the system becomes interference limited. On the contrary, the 
performance under Jakes' model degrades gracefully as the 
user mobility (Doppler bandwidth) increases. 
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